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Collective Excitations of Bose-Einstein Condensates 
in a Double- Well Potential 
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We investigate collective excitations of Bose-Einstein condensates at absolute zero in a 
double-well trap. We solve the Bogoliubov equations with a double-well trap, and show that 
the crossover from the dipole mode to the Josephson plasma mode occurs in the lowest energy 
excitation. It is found that the anomalous tunneling property of low energy excitations is 
crucial to the crossover. 
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1. Introduction 

One of the most intriguing features in Bose-condensed systems is the coherence of the 
macroscopic quantum phase of the order parameter. A Bose-Einstein condensate behaves as 
a quantum mechanical wave with a macroscopic quantum phase. The macroscopic quantum 
phase difference of a condensate trapped in a double-well potential was measured by means of 
the matter wave interference both in destructive^ and nondestructive'^ ways, and the coherence 
of the Bose-Einstein condensate was clearly verified. 

The Josephson effect, which was first predicted for superconductor tunnel junctions,^ also 
clearly exhibits the coherence of a Bose-Einstein condensate. The two mode approximation is 
useful to describe the Josephson-like dynamics of a weakly linked Bose-Einstein condensate 
in a double- well potential (see, e.g., refs. 4-7 and references therein). Using the two mode ap- 
proximation, various types of the Josephson-like dynamics have been predicted, which include 
Josephson plasma oscillation and nonlinear self-trapping.^"^ Recently, the Josephson plasma 
oscillation and the nonlinear self-trapping have been observed by Albiez et al. in a condensate 
of atomic gases trapped in a double-well potential.^ Because of the high tunability of system 
parameters, such a system provides an ideal stage for observing the Josephson-like dynamics. 

Since the Josephson plasma mode is a small amplitude oscillation from a static equi- 
librium, it can be treated as a kind of Bogoliubov excitations.^ When the potential barrier 
separating the condensate does not exist, the lowest energy excitation is the dipole mode. In 
contrast, when the potential barrier is so strong that the two mode approximation can be 
valid, the lowest excitation energy accords with the Josephson plasma energy. In other words, 
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the crossover from the dipole mode to the Josephson plasma mode occurs in the lowest energy 
excitation. 

Salasnich et al. numerically solved the Bogoliubov equations in a double-well trap with 
a harmonic confinement and a Gaussian potential barrier to obtain the collective excitation 
energies. "^'^ They found that the lowest excitation energy decreases as the barrier strength 
increases; however they did not refer to the relation between the reduction of the lowest 
excitation energy and the crossover to the Josephson plasma mode. The excitation energy in 
a double-well potential can be determined by a solution of a scattering problem of Bogoliubov 
excitations through the potential barrier used for preparing the double-well potential. Kagan 
et al. studied the scattering problem of Bogoliubov excitations, and predicted that a potential 
barrier is transparent for low energy excitations, which is called anomalous tunneling}^ It is 
expected that we can relate the anomalous tunneling to the crossover to the Josephson plasma 
mode, when the Josephson plasma mode is regarded as a Bogoliubov excitation. 

In the present paper, we study collective excitations of Bose-Einstein condensates at tem- 
perature T = in a double-well potential, and discuss the crossover from the dipole mode 
to the Josephson plasma mode. Adopting a box-shaped double-well potential, we analytically 
calculate the lowest excitation energy and the Josephson plasma energy for arbitrary values of 
the barrier strength. We show that the lowest excitation energy asymptotically approaches the 
Josephson plasma energy as the barrier strength increases. We also find that the anomalous 
tunneling determines the region of the barrier strength where the crossover occurs. More- 
over, we numerically calculate the lowest excitation energy in case of a double-well potential 
consisting of a harmonic confinement and a Gaussian potential barrier. It is shown that the 
mechanism of the crossover is valid also in this experimentally accessible trap. 

The present paper is organized as follows. In §2, solving the Gross-Pitaevskii equation with 
a box-shaped double-well potential, we analytically calculate the condensate wave function 
and the Josephson plasma energy. In §3, solving the Bogoliubov equations with the box- 
shaped double-well potential, we analytically obtain the excitation spectra of condensates. 
In §4, focusing our attention on the lowest energy excitation, we discuss the crossover from 
the dipole mode to the Josephson plasma mode. In §5, we numerically calculate the lowest 
excitation energy for a condensate in an experimentally accessible double-well potential. In 
§6, we summarize our results. 

2. Condensate Wave Functions and Josephson Plasma Energy in a Box-Shaped 
Double- Well Potential 

We first consider a Bose-Einstein condensate in a box-shaped trap which consists of a 
radial harmonic confinement and end caps in the axial direction (the x axis). We assume 
that the frequency lo± of the radial harmonic potential is large enough compared to the 
excitation energy for the axial direction. Then, one can justify the one-dimensional treatment 
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Fig. 1. Schematic picture of a condensate in a box-shaped double- well potential. Black solid line 
represents the condensate wave function. 

of the problem. Such a configuration was realized in a recent experiment. Setting a potential 
barrier in the center of the trap, a double-well potential is created. Adopting rigid walls and 
a (5-function potential barrier as the end caps and the barrier respectively, the double-well 
potential is written as 



where a is the size of a well. 

Since we are mainly interested in collective excitations of a double- well trapped condensate 
at the absolute zero temperature, our formulation of the problem is based on the Bogoliubov 
theory. The Bogoliubov theory consists of the time-independent Gross-Pitaevskii equation and 
the Bogoliubov equations.^ In this section, we solve the time-independent Gross-Pitaevskii 
equation: 



and obtain the condensate wave function ^'0(2;). Here /u is the chemical potential and m is 
the mass of an atom. Since the radial confinement is harmonic and sufficiently tight, the 
coupling constant is affected by the harmonic oscillator length a± of the radial confinement 
as g = where a, is the s-wave scattering length. ""^^ 

In §4, we will discuss the crossover from the dipole mode to the Josephson plasma mode 
by comparing the lowest excitation energy to the Josephson plasma energy. Using solutions 
of the Gross-Pitaevskii equation, one can calculate the Josephson plasma energy. In order to 
calculate the Josephson plasma energy, we need to obtain not only the lowest energy symmetric 
solution (^'o™) but also antisymmetric (^'g'^) solution of eq. (2). Assuming that the system 
size is much larger than the healing length ^ = —fi=, the condensate wave function near the 



center of each well is not affected by the potential barrier and the rigid wall. Then, one can 
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approximately obtain the symmetric and antisymmetric solutions of eq. (2) as 



tanh ( ^ 



9 \ tanh (^) , I < |x| < a, 



*o» = a/^x<; 



tanh ( ^ ) , kl < §, 



sgn(a;) tanh , § < |x| < a, 

where /Xsm(an) Csm(an) is the chemical potential and the healing length of the symmetric 
(antisymmetric) state. A constant xq reflects value of the condensate wave function at x = 0, 
which is determined by the boundary conditions 



*o(-0) = ^'oC+O), (5) 



dx 

as 



dx 



x=+0 "' 



'"""VolO). (6) 



ssm -^Msmssm 

We next calculate the Josephson plasma energy ejp = EjEq , which is easily derived 
from the Josephson Hamiltonian: 

Hj = - Ej cosif, (8) 

where k and if represent the population difference and the phase difference between the two 
wells. ^ The Josephson coupling energy Ej expresses the overlapping integral of condensate 
wave functions in the two wells, and the capacitive energy Eq is proportional to inverse of 
the compressibility. They are defined as 

El = (9) 

Ea = 4^, (10) 

where E^^f^g^^-^ is the mean field energy of the symmetric (antisymmetric) state.^ Hence, we 
need to calculate the chemical potential and the mean field energy in order to obtain the 
Josephson plasma energy. 

The chemical potential is related to the number Nq of condensate atoms by the normal- 
ization condition: 

dx\^o{x)f = No. (11) 

Substituting eqs. (3) and (4) into the normalization condition, we derive the relations between 
the chemical potentials and Nq: 

/^sm A 2^ + ^!f£tanh^^ = 1, (12) 
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Fig. 2. Chemical potentials /ism (solid line) and /ian (dashed line) as functions of the potential 
strength Vq are shown. 
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where ng = ^ is the averaged density of the condensate. One can obtain approximate so- 
lutions of eq. (12) in the limits of Vq » gn^^Q and Vq ^ gno£,o, where = ^mgno ■ When 



Vq ^ gno^o, we expand eq. (12) into power series of ^ and ^"y^° , and obtain 



/fsm _ ^ 2go ^ 2gg 
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In a similar way, when Vq <C gn^^Q^ we expand eq. (12) into power series of ^ and — ^ 



9'^0?0 



obtain 
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Meanwhile, the solution of eq. (13) is 



(14) 
and 

(15) 
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+ 



+ 



. . (16) 
In eqs. (14), (15) and (16), we express the expansions up to the third order of the small 
parameters. In Fig. 2, we show the chemical potentials ^sm and /Xan as functions of the potential 
strength Vq. While fia_n is constant, fism increases monotonically and approaches the value of 
/Xan as Vq increases. This means that the symmetric state and the antisymmetric state are 
degenerate in the strong potential limit Vq ^ oo. 

The mean field energy E can be calculated using the expression 



E 



dx^l{x) 



(17) 



5/14 



J. Phys. Soc. Jpn. 



Full Paper 




20 



40 60 80 

Fo/(gno4) 



100 



Fig. 3. Josephson plasma energy ejp as a function of the potential strength Vq is shown. 



In the limit of Vq ^ 9^0^,0, one can approximately obtain 
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(18) 
(19) 



where -Esj^(an) is the mean field energy of the symmetric (antisymmetric) state. 

Substituting the chemical potential of eq. (14) and the mean field energies of eqs. (18) and 
(19) into eqs. (9) and (10), one obtains 



griQ 

Eq 
gno 



4Voa 



1 + — 

a 



when Vq ^ gn^^Q. These equations yield the Josephson plasma energy 



£JP 

gm 



1 + 



(20) 
(21) 

(22) 
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The Josephson plasma energy is shown in Fig. 3, as a function of the potential strength. In 
§4, we will see that the lowest excitation energy coincides with this Josephson plasma energy 
for the sufficiently strong potential. 

3. Excitation Spectrum in a Box-Shaped Double- Well Potential 

The aim of this section is to calculate the excitation spectra of condensates in a box-shaped 
double-well potential. The excitations of a condensate correspond to the small fiuctuations of 
the condensate wave function ^'o(x,t) around the stationary solution of eq. (2), 



5^'o(x,t) = e 'ft* ^ ( u,j{x)e 



v*{x)e 



(23) 
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where £j is the energy of the excitation in an eigenstate labeled by j. The wave functions 
{uj{x),Vj{x)Y of the excitations fulfill the Bogoliubov equations,^ 



Ho -g^oixf 




(24) 



^' = ~2^^d^~^^^^^^^ 25|^'o(x)|^ (25) 
Solving these equations with the condensate wave function of cq. (3), we shall obtain the 
excitation spectrum. Hereafter the chemical potential n denotes the chemical potential //gm of 
the symmetric state. 

In order to relate the problem to tunneling properties of the Bogoliubov excitation, we 
separate the system into three regions, or —a < x < — |, |x| < | and | < x < a. We solve 
the Bogoliubov equations analytically in each region, and derive the equation to determine 
the excitation spectrum by connecting each solution smoothly. 

At first, we find solutions of the Bogoliubov equations in the central region (|x| < |). 
There exist two independent solutions with a same excitation energy e, corresponding to 
two types of scattering process. One solution tp^x) = (u'(x), f'(x))* describes the process 
where a Bogoliubov excitation comes from left, and the other solution ip'^{x) = (u'^ {x) , v'^ {x))^ 
describes the process where a Bogoliubov excitation comes from right. The solution tp^x) is 
written as 




-p-X 




(26) 



U2 ' 
V2 J 

Va J 

where p± = \j'^{^\fi^ + ^'^ =f70 satisfies the Bogoliubov spectrum for a uniform sys- 
tem.^^'^^'^^ The coefficients r, 6, and c are the amplitudes of reflected, left localized, trans- 
mitted, and right localized components, respectively. All the coefficients can be determined 
by the boundary conditions at x = 0. Thus, we obtained ip^x), and we can also find ^''(x) in 
the same way. 

Expanding \t\ and the phase S of t around e = 0,^^ one can analytically obtain 

\t\ - i-<»(j)'. (27) 

S ~ (28) 
A* 

The coefficients a and (3 are 

2{Vo-,,0{Vi+uVi +2u{fiO ^ - 4(/xO=^) + 9(m^Fo) 



^ ^ V^ + ,yV0-3fi^l^ + 6if,O^ 

2/x^z^ 



7/14 



J. Phys. Soc. Jpn. 

where 



Full Paper 



^ = ^Jvs+m?. (31) 

It is obvious from eqs. (27) and (28) that the transmission coefficient approaches unity 
and the phase shift of t approaches zero as the energy is reduced to zero; Kagan et al. called 
this behavior anomalous tunneling}^ This expression is adequate for arbitrary values of the 
potential strength Vq. Furthermore, when £ <C /i and Vq S> fx^, one can approximately obtain 
the reflection and transmission amplitudes/^ 



e^Vo 



r = 



eVo - £/x^ + i 
eVo — en^ + iii^^ ' 



(32) 



(33) 



In eq. (33), we can see that the peak around e = in |tp has the Lorentzian shape with 
half width Ae ~ Such an anomalous tunneling property of the Bogoliubov excitations 
essentially affects the crossover from the dipole mode to the Josephson plasma mode of the 
lowest energy excitation, as discussed in the next section. 

We can write a general solution of the Bogoliubov equations in the central region as a 
linear combination of ^'(x) and ^''(x): 

(34) 

One can obtain an analytical solution of the Bogoliubov equations also in the left and 
right side regions (| < |x| < a). It is 



t/j{x) = r]'^\x) +C,ilf{x), \x\ < ^. 



w 



side 



(X) = 



v^'^'^'^ix) 




+ 1 - tanh^ 



\x\ 



2e 



2£ I 2 



tanh 



tanh 



X — a 



cosp+(|x| — a) 

sinj3_|_(|x| — a) 
cos;j+(|x| — a) 



(35) 



sinp+(|x| — a) 



(36) 



Using the solutions of eqs. (34), (35) and (36), one can construct solutions in all regions: 



■0(x) = < 



^^side(-^) -a < X < -§, 

rj->p'{x) + C'i/'''(x) |x| < §, 



(37) 



GV'''^^(x) § < X < a, 

where 1] = ( and F = G are satisfied for collective excitations with even parity ot r) = —( 
and F = —G for those with odd parity. By imposing the boundary condition at |x| = |, the 
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Fig. 4. Excitation spectra in the double-well traps with Vb = 0, Vb = 10 and Vb = 100 are shown, 
where a = lOfo- 

solutions can be smoothly connected. As a result, one obtains the equation to determine the 
excitation spectrum: 

(r ± t)exp (i(2p+a + 7)) = 1, (38) 

where 

7 = tan-i ( ] . (39) 

The positive (negative) sign in the left-hand side of eq. (38) represents the excitations with 
even (odd) parity. One can see from eq. (38) that tunneling properties of the Bogoliubov 
excitations through the potential barrier affect the excitation spectrum explicitly. Since the 
reflection and transmission amplitudes satisfy relations 

+ |r-|2 = 1, (40) 

and 

t = |t|e^^ r = ±i\r\e'^, (41) 

one obtains 

|r ± t| = 1. (42) 
Using eq. (42), one can rewrite eq. (38) in terms of the phase of the left-hand side as 

2p+a + 7 + (/>± = 27rn, (43) 

where (j)± is the phase of r ±t. 

Solving eq. (43) for excitation energy e, we obtain the excitation spectra which consist 
of discrete energy levels as shown in Fig. 4. Parity of quantum numbers expresses parity of 
the wave function of the collective excitations. In the strong potential limit Vq ^ 00, the 
excitations labeled by 21 and 21 + 1 are degenerate, because the condensate is completely 
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divided into each well. Hence, as the potential barrier becomes stronger, energy of the even 
parity excitation labeled by 21 increases and energy of the odd parity excitation labeled by 
21 + 1 decreases so that both energies approach. The change of energies of the odd parity 
excitations is more pronounced than that of the even parity excitations, and such a tendency 
is qualitatively consistent with the numerical results in the case of the double-well traps with 
harmonic confinement.^'' Moreover, we see that the less the excitation energy is, the slower 
the change according to growth of the potential barrier is. This implies that the anomalous 
tunneling behavior of low energy excitations buffers the effect of the potential barrier. 

4. Crossover from Dipole Mode to Josephson Plasma Mode 

In this section, we shall discuss the effect of the anomalous tunneling on the lowest energy 
excitation. Solving eq. (43) with n = and odd parity, we obtain the lowest excitation energy 
as a function of the potential strength, which is represented by the solid line in Fig. 5. In 
order to elucidate the effect of the anomalous tunneling on the lowest excitation energy eiow, 
we calculate it analytically in two limits. 

On the one hand, when S> ^, substituting the low energy expansion of the trans- 

mission amplitude of eqs. (27) and (28), one obtains 



gno 2a 



(44) 



a \2 2 V 2 4^710^0 
which is plotted in Fig. 5 by the dotted line. Obviously, the influence of the potential barrier 
on the lowest excitation is relatively small, because the correction of £iow due to the potential 
barrier is included in the small term of O(^) in eq. (44). On the other hand, when <C ^, 
one can obtain an approximate expression of the lowest excitation energy from eqs. (32) and 
(33), 



gno V V 2a / 

This energy coincides with the Josephson plasma energy of eq. (22), and this result explicitly 
justifies treatment of the Josephson plasma oscillation by the two mode approximation for a 
sufficiently strong potential barrier. In Fig. 5, one can see that the crossover from the dipole 
mode to the Josephson plasma mode occurs. 

It is crucial to understand the physics of the crossover that two energy scales of, ^ and 
^^"{/p ^° ) determine whether the lowest energy excitation is the Josephson plasma mode or the 
dipole mode. The former is comparable to the dipole mode energy. The latter is comparable to 
the half- width Ae of the peak of the transmission coefficient In other words, the potential 
barrier is almost transparent for the excitations with energy less than ^^"y^ ^° • The crossover 
is determined by whether the anomalous tunneling is effective or not for the lowest energy 
excitation. When ^ <C , the potential barrier is almost transparent for the lowest energy 
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Fig. 5. Black solid line is the lowest excitation energy for a — lO^o as a function of Vq. Dotted line is 
the corrected dipole mode energy of eq. (44), while the dashed line is the Josephson plasma energy 
shown in Fig. 3. Gray solid line is the half-width of the peak of the transmission coefHcient. 

excitation due to the anomalous tunneling; . the dipole mode is hardly affected by the poten- 
tial barrier. As the potential strength Vq increases, contribution of the anomalous tunneling 
diminishes gradually and the lowest energy excitation comes to be significantly affected by 
the potential barrier Consequently, when ^ ^ ^t^i the lowest energy excitation becomes 
the Josephson plasma mode. Figure 5, where the half- width of the anomalous tunneling (gray 
line) divides the lowest excitation energy into regions of the dipole mode and the Josephson 
plasma mode, clearly confirms such an interpretation of the crossover. 

5. Realistic Double- Well Potential 

In the preceding sections, we considered the case in the double-well potential of eq. (1) 
with rigid walls and a (5-function potential barrier. In this section, we numerically solve the 
Gross-Pitaevskii equation and the Bogoliubov equations in case of a realistic double-well 
potential. 

Considering a double-well potential consisting of the magnetic trap and the blue-detuned 
laser beam, which is realized in the experiment of ref. 1, we adopt a double- well potential with 
a harmonic confinement and a Gaussian potential barrier: 



where Ux is the frequency of the harmonic potential. The height Uq and the width a of the 
potential barrier can be controlled in experiments by varying the intensity and the aperture 
of the laser beam, respectively. 

With use of the Bogoliubov theory, one can discuss the crossover to the Josephson plasma 
mode in the same way as in the case of the box-shaped double-well potential. We can numeri- 
cally find the lowest symmetric and antisymmetric solutions of the Gross-Pitaevskii equation 
with the double- well potential of eq. (46) ; accordingly we can calculate the Josephson coupling 




(46) 
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Fig. 6. Closed circles represent the lowest excitation energy, and open circles represent the Josephson 
plasma energy. Triangles represent the half- width of the peak of the transmission coefficient. We 
consider a condensate of ^'^Na atoms. The values of parameters are as follows; the number of the 
condensed atoms is A^o = 3000, the frequency of the radial confinement is uj± — 250 x 2tt Hz, the 
frequency of the axial confinement is lu±_ = 10 x 27r Hz, the s-wave scattering length is — 3nm 
and the width of the potential barrier is = 1.33 /im. 



energy Ej and the capacitive energy Eq using the expressions of eq. (9) and eq. (10). As a 
result, we can obtain the Josephson plasma energy ejp. Furthermore, substituting the lowest 
symmetric solution of the Gross-Pitaevskii equation into the Bogoliubov equations, we can 
numerically calculate the lowest excitation energy. 

Open circles and closed circles in Fig. 6 represent the lowest excitation energy and the 
Josephson plasma energy, respectively. In Fig. 6 one can see that the lowest excitation energy 
reduces as height of the potential barrier increases; finally it accords with the Josephson 
plasma energy for the sufficiently high barrier. Thus, the crossover from the dipole mode to 
the Josephson plasma mode is confirmed in the case of a realistic double-well potential. 

In order to elucidate the role of the anomalous tunneling in the crossover, we next estimate 
half-width Ae of the peak of the anomalous tunneling for a Gaussian potential barrier. In case 
of a rectangular potential barrier, Kagan et al. have approximately obtained an analytical 
expression for the transmission coefficient.^^ Using the expression, one can calculate the half- 
width. 



(47) 



where 



1 /"a 2m 

'^^^ H I d^^y (^barrier (a;) - ^). (48) 

In the region \x\ < ^, the potential barrier Vbarrier(2;) is larger than the chemical potential; 
this means that the points |x| = ^ are the classical turning points. We estimate the half-width 
Ae for a Gaussian potential barrier by means of the expression of eqs. (47) and (48). 
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Triangles in Fig. 6 represent the half-width Ae of the peak of the transmission coefficient. 
When the dipolc mode energy is much smaller than the half-width, namely fkVx <C Ae, 
the anomalous tunneling is effective for the lowest excitation energy. Accordingly, the lowest 
excitation is hardly affected by the potential barrier. When fkOx ^ the lowest excitation 
energy remarkably changes because the excitation begins to perceive the potential barrier. 
When fujx ^ Ae, the anomalous tunneling is no longer effective; consequently, the lowest 
excitation becomes the Josephson plasma mode. Therefore, the interpretation of the crossover 
discussed in the previous section is also appropriate for the double- well potential of eq. (46). 

6. Conclusions 

In summary, we have studied collective excitations of a condensate in a double-well po- 
tential. We have analytically solved the Bogoliubov equations with a box-shaped double-well 
potential, and analyzed the excitation spectrum. In the lowest excitation, it has been found 
that the crossover from the dipole mode to the Josephson plasma mode occurs as the po- 
tential barrier separating the condensate becomes strong. The crossover is dominated by the 
anomalous tunneling behavior of the Bogoliubov excitations. We have numerically calculated 
the lowest excitation energy and the Josephson plasma energy for a condensate in a realistic 
double-well potential. We have confirmed by numerical calculation that the mechanism of the 
crossover is valid also in the case of more realistic potential. 

While only condensates at T = has been considered in our calculations, the Joseph- 
son plasma oscillation of condensates in finite temperature is known to exhibit dissipative 
behaviors due to the presence of the thermal depletion.^ It will be interesting to study the ex- 
citations of the condensates at finite temperature in the double-well potential, and to discuss 
the relation between the damping of the collective oscillation and the anomalous tunneling of 
the Bogoliubov excitations. 
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